Study on the material source and enrichment mechanism of REE-rich phosphorite in Zhijin, Guizhou

Rare earth element (REE)-rich phosphorite in the Guizhou region mainly exists in the Doushantuo Formation and Gezhongwu Formation in early Cambrian strata, which are some of the important strata containing phosphorite resources in China. The early Cambrian Zhijin phosphorite in Guizhou Province, China, has high rare earth element and yttrium (REY) contents of up to 2500 ppm, with heavy REY (HREY) contents accounting for ~ 30% of the total REY contents. However, the specific controlling source and environment of phosphorite (especially the REEs in Zhijin phosphorite) are still unsolved. Through field geological investigations; mineralogical, geochemical, Sr–Nd isotope analyses; and tectonic characteristics, the material source, sedimentary environment and seawater dynamics of REEs in phosphorite are studied. It is considered that the REEs enriched in the Zhijin phosphorite are mainly affected by precipitation from hydrothermal fluid. Moreover, from the late Ediacaran to the early Cambrian, the depositional environment from the bottom to the top of the water tended to be hypoxic, and the activity of hot water fluid tended to be strong. The change in redox conditions is closely related to the rise and fall of sea level. Combined with the tectonic background, these results show that the weakly oxidized environment may be an important factor controlling the enrichment of REEs. The enrichment of REEs may be closely related to volcanic hydrothermal activity, later diagenesis and seawater dynamics.


Background
Regionally, Guizhou Province is part of the Yangtze plate, which is mainly composed of phosphorite from the Ediacaran Doushantuo Formation and Cambrian Gezhongwu Formation.The large phosphorite deposits are mainly distributed in Weng'an, Fuquan, Kaiyang, and Zhijin.In addition, medium-sized phosphorite deposits are distributed in Songtao, Jinsha, Duyun, Zunyi and other areas (Fig. 1).According to the new regional geological map of Guizhou Province, the typical phosphorite of the Cambrian Gezhongwu Formation is located in southeastern Zhijin County, Guizhou Province (Fig. 2).The typical phosphorite of the Ediacaran Doushantuo Formation is located in Weng'an, Guizhou Province (Fig. 3).
The tectonic lines (including faults and folds) in the Zhijin phosphorite deposit are mainly spread from north to east, and the tectonic deformation is mainly dominated by positive faults, reverse faults, and strike-slip faults (Fig. 2).Among them, the Gezhongwu phosphorite area is controlled by the northeast Guohua-Gezhongwu fault, which makes the Cambrian formation spread out along the fault, consistent with the fault and anticline direction.The Damachang phosphorite area is controlled by the Guiguo-Damachang fault.The Maobaichong phosphorite area is controlled by a series of northeast faults.The strata in the Zhijin phosphorite deposit are Ediacaran (Z), Cambrian (ε), Devonian (D), Carboniferous (C), Permian (P), and Triassic (T) (Fig. 2).Among them, the Cambrian Gezhongwu Formation is the formation bearing REE-rich phosphorite.The lower ore layer is composed of dolomaceous phosphorite and phospho-dolomite, which are rich in microorganisms.The upper ore layer is phospho-claystone and siltstone.The Gezhongwu Formation unconformably contacts the underlying Dengying Formation, and the overlying strata are carbonaceous shales of the Niutitang Formation 19 .
The Weng'an phosphorite deposit in Guizhou is located in the southeastern margin of the Yangtze Craton.The phosphorite originates from the nucleus and wing of the Baiyan-Gaoping anticline.The anticline saddle part is truncated by the north-south Xiaoba fault, with the Baiyan anticline in the north and the Gaoping anticline in the south.Both the anticlines and the Weng'an phosphorite deposit are controlled by the Xiaoba fracture.The strata in the Weng'an phosphorite are Ediacaran (Z), Cambrian (ε), Permian (P), Triassic (T) and Jurassic (J) (Fig. 3).Among them, the late Ediacaran Doushantuo Formation is the formation bearing phosphorite.It is in false contact with the underlying Nantuo Formation and conforms to the overlying Dengying Formation.The lower ore layer is sandstone phosphorite, and the upper ore layer is dolomaceous phosphorite 19 .

Sample collection
Fourteen samples were collected from the Gezhongwu Formation of the Zhijin phosphorite deposit (Fig. 4a-c).For comparison, three additional samples were collected from the Doushantuo Formation of the Weng'an phosphorite deposit.Sample features are shown in Table 1.

Petrographic analysis
The structures in phosphorites are banding, lamination, massive, brecciated, spherical, cataclastic, clumpy and earthy.Among them, the structures in phosphorites in the Doushantuo Formation are mainly banding (Fig. 4d), massive (Fig. 4e), brecciated (Fig. 4f), and earthy (Fig. 4g).Banding is the most common structure in phosphorite, with a regular dark gray and light white pattern and a high ore grade.The massive structure is brownish yellow, dense and hard.The mineral particles are evenly distributed, undeveloped or rarely intermittent in dolomite layers.The brecciated structure is dark gray breccia cemented by light gray dolomite, with particle sizes of 2-3 mm, indicating that the primary phosphorite was broken and cemented by dolomite after migration.The earthy structure is earth-yellow, fragile and soft.This structure is always found in weathered phosphorite.
The structures in phosphorites in the Gezhongwu Formation are mainly laminations (Fig. 4h), followed by clumpy (Fig. 4i), cataclastic (Fig. 4j,k) and spherical structures (Fig. 4l).The laminations are intermittently developed with flocculent, brecciated and clumped structures interspersed in dolomite, and a few layers have thicknesses of less than 1 mm.The clumpy structure is an ellipsoid with particle sizes of 3-5 cm.The clumps are pure phosphorite and quartz.The cataclastic structure shows irregular breccia with particle sizes of 2-5 mm and is filled by dolomite veins.The spherical structure is spherical, with gravel and sandstone phosphorite cemented by dolomite, and the particle sizes of sand are 3-5 mm.
Parallel bedding is present in the Doushantuo Formation phosphorites of Weng'an (Fig. 4m), while parallel bedding and cross-bedding are present in the Gezhongwu Formation phosphorites of the Zhijin deposit (Fig. 4n).Wave bedding is present in some areas of the Gezhongwu Formation phosphorites (Fig. 4o).
This shows that the Zhijin phosphorites were subjected to stronger seawater dynamic action in the process of sedimentation.The particle size of the Zhijin Gezhongwu Formation phosphorites is larger than that of the Weng 'an Doushantuo Formation, and the structural maturity of the Zhijin Gezhongwu Formation phosphorites is higher.Compared with Weng'an phosphorites, Zhijin phosphorites were subjected to greater mechanical power from seawater, and the seawater depth was shallower.
Microscopic rock and mineral identification was completed in the Earth Science Laboratory of Zhejiang University, with a magnification of 5-50 times.The mineral composition and texture of the rocks were observed under single polarized light, orthogonally polarized light and reflected light.Phosphorite from the Zhijin and Weng'an deposits consists of phosphate, dolomite, feldspar, quartz, calcite, mica and other minerals.The mineral texture mainly includes internal clastic and bioclastic structures and amorphous collophanite.
Mineralogical characteristics show that Weng'an phosphorite is common in internal clastic structures, mainly spherical granular (Fig. 5a), fine-grained (Fig. 5b,c), granular (Fig. 5d), irregular block (Fig. 5e), and more terrigenous clastic material.Among them, the most common structure is the spherical granular structure.It has large particle sizes of approximately 200-500 μm and is ellipsoid, spherical, spindle, highly rounded and has an obvious boundary with dolomite.The fine-grained structure shows that the phosphate is filled with smaller particles of quartz or dolomite internally, which formed by absorbing the quartz and dolomite debris that formed in seawater during the phosphorite precipitation process.Some of the shells are wrapped around the dolomite and are approximately 10-20 μm thick, forming a granular structure.In addition, internal debris also appears in Figure 1.Tectonic sketch map and phosphorite distributions in Guizhou 19 (modified from reference 19 and created it with Adobe Illustrator CS6.http:// www.adobe.com/ produ cts/ illus trator.html).
the irregular block structure.This shows that after phosphorite precipitation, the phosphorite was not damaged by strong seawater action or weaker seawater and was directly buried for mineralization.
However, the Zhijin phosphorite is common in bioclastic structures, such as conodonts (Fig. 5f), zones (Fig. 5g), and short columns (Fig. 5h), and a small amount of amorphous collophanite is also found (Fig. 5i).The conodont structure shows an open fan shape on one side and a closed cone on the other side, with particle sizes of 100-300 μm.Most conodont biological debris appears porous and relatively loose.The zonal structure shows ribbon-shaped biological debris, particle sizes of 10-1000 μm, and relatively large length and width ratios, mostly 8:1-10:1, with the characteristics of bending.Most zonal structures have rough surfaces, and some fragments have dense phosphorite inside.Short columnar structures are mainly short, cylindrical and partially dense.The particle sizes of columnar biological debris are 50-100 μm, and most of the debris are concave and uneven.The amorphous collophanite is in an unfixed form, cemented with dolomite, or distributed between dolomite crevices.
Mineralogical analysis indicated that during the Ediacaran period, there was no biological activity involved.After the deposition of phosphorite, the geochemical properties of pore water were stable, so the spherical phosphorite fragments of the Doushantuo Formation were preserved intact.During the early Cambrian period, algae and microorganisms proliferated and degraded after death.The released P formed phosphorite in pore water, resulting in the redistribution and differentiation of REEs and P. Therefore, some phosphorites in the Gezhongwu Formation exhibit destroyed biological residues, with amorphous collophanite and dolomite cemented together.The microorganisms were quickly buried after death, also leading to phosphorylation at the sedimentary interface.

Chemical analysis
The 17 fresh samples were selected and crushed in a ceramic jaw crusher and powdered to ~ 200 mesh in an agate mill for whole-rock elemental and isotopic analyses.(modified from reference 20 and created it with Adobe Illustrator CS6.http:// www.adobe.com/ produ cts/ illus trator.html).
Vol.:(0123456789) www.nature.com/scientificreports/Major element oxides were analyzed by the X-ray fluorescence (XRF) method (Rigaku RIX2100) using Liborate glass disks.BCR-2 and GBW07105 were external standards, and the analytical accuracy and precision were both better than 5%.
Trace element concentrations were determined using an Agilent 7500a inductively coupled plasma mass spectrometry (ICP-MS) instrument.All samples were dissolved in a sealed high-temperature and high-pressure bomb using an equal proportion of superpure HF and HNO3.Details on the analytical procedure can be found in Liu et al. 22,23 BHVO-2, AGV-2, BCR-2 and GSP-1 were used as standards, and the analytical accuracy was mostly better than 5% (relative error).
The REEs of all samples were measured via inductively coupled plasma atomic emission spectrometry (ICP-AES) (Thermo iCAP-6300).The mineral composition and content were analyzed by FESEM-EDS (FESEM; ZEISS Sigma-500, EDS; BRUKER XFlash 6|60) with an advanced mineral identification and characterization system (AMICS).The experimental conditions were as follows: the acceleration voltage was 20 kV, the resolution was 0.8 nm, and the probe current was 40-100 nA.The above analyses and testing were performed at the Hongchuang Testing Center of Nanjing University.δCe and δEu were calculated by δCe = Ce N /(La N × Pr N ) 1/2 and δEu = Eu N /(Sm N × Gd N ) 1/2 , respectively, and N is the average of the post-Archean Australian shale standard value. 24All the original values for the major elements, trace elements and rare earth elements of the 17 samples are shown in Table 2.
The Sr-Nd isotope separation test of rock samples was completed on the MC-ICP-MS from the laboratory of Nanjing Hongchuang Detection Technology Co., Ltd.The Rb-Sr and Sm-Nd isotope separations of wholerock samples included the processes of sample preparation, chemical separation, and testing (Table 3).The rock powder was dried in an oven and then weighed.Similar to the method of treating samples by ICP-MS, the acid Figure 3.The tectonic geological picture of Weng'an phosphorite deposit 21 (modified from reference 21 and created it with Adobe Illustrator CS6.http:// www.adobe.com/ produ cts/ illus trator.html).dissolution method was used to dissolve the sample.Rb-Sr and Sm-Nd isotope separation was carried out using the conventional cation exchange resin method.The type of resin used for Rb-Sr separation was AG-SOW-X8, and HDEHP resin was used for Sm-Nd separation.For detailed separation and test procedures, see Liu et al. 25

Results
The statistics of the major elements, trace elements, rare earth elements and different element ratios are calculated and listed in Table 4.It is noteworthy that the major element contents in the Doushantuo Formation vary greatly from those in the Gezhongwu Formation of the Zhijin deposit.The phosphorite from the Doushantuo Formation of the Weng'an deposit is characterized by a high CaO content (the average value is 52.93%), low SiO 2 content (average value of 1.49%), Al 2 O 3 content of 0.32%, TFe 2 O 3 content of 0.25% and (K 2 O + Na 2 O) content of 0.43%.Moreover, the content of MgO is low, with an average value of 1.31%.However, the content of P 2 O 5 (the average value is 36.83%) in the Doushantuo Formation of the Weng'an deposit (the average value is 16.08%) is much higher than that from the Gezhongwu Formation of the Zhijin deposit.

Characteristics of trace elements
The PAAS normalization of trace elements (PAAS data are from Taylor and Mclennan 26 ) for the samples collected in this study indicated that the enrichment degree of trace elements differs in different strata from the Zhijin and Weng'an phosphorite deposits.Among them, the dolomite from the Dengying Formation is rich in Sr, Ba and U; phosphorites from the Gezongwu Formation of the Zhijin deposit are characterized by enrichment in Zn, Sr, and U; and carbonaceous shales from the Niutitang Formation are enriched in V, Ni, Zn, Sr, and U, as shown in Fig. 6.The enrichment degrees of Sr and Ba vary in different strata in the two deposits.Moreover, the enrichment degrees of different trace elements in the Gezhongwu Formation vary partially.www.nature.com/scientificreports/

Geochemical characteristics of rare earth elements
REEs in phosphorite deposits from the Zhijin and Weng'an deposits in Guizhou vary greatly in different strata (Fig. 7).The Gezhongwu Formation contains the highest REE content in this study, with an average ΣREE content of 801.57ppm and values varying from 508.25 to 1409.6 ppm.LREE/HREE ratios vary from 0.81 to 2.05, with an average of 1.33.The PAAS normalization pattern is left-leaning, with cap-type distribution characteristics and a significant negative Ce anomaly (varying from 0.03 to 0.52, with an average value of 0.22).δEu values vary from 1.03 to 1.44, with ZJ-5 exhibiting a high δEu of 1.44.Most samples have δEu values of approximately 1, exhibiting a positive Eu anomaly.
Dolomite from the Dengying Formation contains highly variable ΣREEs from 5.99 to 527.19 ppm, with an average value of 236.29 ppm.The LREE/HREE ratios range from 0.78 to 4.92, and the average value is 1.92, exhibiting enrichment in LREEs.The samples have negative Ce anomalies (varying from 0.05 to 0.99, with an average of 0.41); in addition, most samples exhibit positive Eu anomalies, except sample ZJ-4.
The carbonaceous shale from the Niutitang Formation contains relatively low ΣREEs, varying from 114.22 to 259.61 ppm, with an average content of 196.65 ppm.LREE/HREE ratios show high fractionation between LREEs and HREEs, varying from 1.92 to 5.96, with an average of 3.74, indicating distinct LREE enrichment.It is noteworthy that the shale shows significant positive Ce anomalies and slight negative Eu anomalies.
Comparatively, the samples from the Gezhongwu Formation of the Weng'an deposit contain very low ΣREE contents varying from 44.38 to 177.30 ppm, with an average value of 136.32 ppm.The LREE/HREE ratios range from 1.34 to 2.25, with an average ratio of 1.83, indicating slight enrichment of LREEs.
Most samples in this study exhibit similar PAAS normalization patterns, except for samples from the carbonaceous shale of the Niutitang Formation.

Sr-Nd isotopic composition
The values of the whole-rock Sr-Nd isotope test data for the samples are shown in Table 3.The results of the analysis showed that the ( 87 Sr/ 86 Sr) i values of the Doushantuo Formation range from 0.709997 to 0.710200, and the ε Nd (t) values range from -7.67 to -6.99.The ( 87 Sr/ 86 Sr) i values of the Dengying Formation range from 0.709944 to 0.710165, and the ε Nd In Fig. 8a, the samples plot between the upper crust and the mantle area, showing the characteristics of mantle mixing.The average TDM2 value of the Doushantuo Formation is 2861 Ma, while the average ε Nd (t) value (− 11.60 to − 10.26) is − 11.09, reflecting that the rock in the source region is associated with partial melting of Neoarchean crustal material.The average TDM2 value (2228-2729) of the Gezhongwu Formation is 2414 Ma, while the average ε Nd (t) value (− 8.85 to − 5.85) is − 7.39, which is higher than that of the samples in the Doushantuo Formation.In Fig. 8a, the samples plot in the ancient lithospheric mantle region, indicating that the source rock may have formed in the process of ancient crustal reconstruction caused by volcanic activity.The average ε Nd (t) values of the Dengying Formation and Niuhootang Formation are − 9.62 and − 8.44, respectively, and the average TDM2 values are 2686 Ma and 2542 Ma, respectively, reflecting that the Dengying Formation contains more ancient crustal material in the source rock.
In addition, in Fig. 8b, a simple simulation calculation of Sr-Nd isotopes indicates that the contribution from ancient crustal material is approximately 30-60%, which is within the range of the red dotted line.The order of the contribution of ancient crust material is the Doushantuo Formation > Dengying Formation > Niutitang Formation > Gezhongwu Formation.Among them, the contribution of the ancient crustal material in the Doushantuo Formation is the highest, accounting for more than 50%.

Material sources analysis
In this study, it is indicated that although the grade of phosphorite from the Doushantuo Formation of the Weng'an deposit is much higher than that from the Gezhongwu Formation of the Zhijin deposit, the REE content in the latter is much higher than that in the former.It is then deduced that the supply of REEs is the key factor controlling the REE content in phosphorite.The identification of REE material sources is thus very important to understand the enrichment mechanism of REEs in sedimentary phosphorite.
Usually, marine deposition is composed of extrinsic material, including terrigenous material carried by river runoff, aeolian or glacial debris, and endogenous materials, including volcanogenic material and related alteration products and hydrothermal materials.In addition, precipitation from seawater is another important material provider.To date, most studies have focused more on the material sources of phosphorus in phosphorite from Southwest China with good ore quality, few studies have been performed on the REE material source, and no consensus has been reached for the REEs enriched in phosphorite from the Gezhongwu Formation.Possible REE material sources include hydrothermal material and related alteration products 6 .Alternatively, at appropriate temperatures, REEs and phosphorite exchange substances in pore water, resulting in enrichment of REEs in phosphorite 30 .Wu et al. believed that the enrichment of REEs may be affected by later diagenetic weathering.The sources of REEs are diverse, and various factors need to be studied in more detail 1 .
However, concerning the uneven distribution of REEs from phosphorites from the Gezhongwu Formation and Doushantuo Formation, it is deduced to be impossible for local seawater to be a main REE supplier with great variation in chemical composition in the two periods.In addition, previous studies reveal that the Kunyang phosphorite, which formed in the same geological period as the Zhijin phosphorite, contains very low REE contents 9 ..The ancient lithospheric mantle data are from Ma et al. 29 .
Vol:.( 1234567890 Fe, Mn, Al and Ti are further used to carry out the source study.It is generally believed that the Al-rich marine sediments are related to the intervention of terrigenous materials, while the enrichment of Fe and Mn is mainly affected by the participation of hydrothermal fluid.Therefore, Al/(Al + Fe + Mn) is often used to evaluate the hydrothermal influence 31 .In this sense, marine sediments related to terrigenous material show Fe/ Ti > 20, (Fe + Mn)/Ti > 25, and Al/(Al + Fe + Mn) > 0.5, and the source of hydrothermal action is indicated by Al/(Al + Fe + Mn) < 0.35 31 .Therefore, Fe/Ti-Al/(Al + Fe + Mn) and (Fe + Mn)/Ti-Al/(Al + Fe + Mn) are used to distinguish the rock provenance in this study.The results show that all samples of the Gezhongwu Formation and two of the samples of the Dengying Formation dolomite plot in the hydrothermal sedimentary area (Fig. 9), which is consistent with previous research conclusions 6 .

Evidence from trace elements
In general, it is believed that Cu, Ni and Zn are primary hydrothermal sources that are also involved in biological productivity 32 .The Co/Zn ratio of hydrothermal-sourced sediments is relatively low, with an average ratio of 0.15, while it is much higher in ferromanganese crusts, at approximately 2.5 33 .In this study, Co/Zn in all samples is very low.The Co/Zn ratios in the Weng'an samples range from 0.01 to 0.02, while the Co/Zn ratios in the Zhijin samples range from 0.01 to 0.20.With the usage of the Zn-Ni-Co ternary diagram (Fig. 10), it is noteworthy that the phosphorite from the Gezhongwu Formation and dolomite from the Dengying Formation are located in the hydrothermal sedimentation area and vicinity, showing distinct involvement of hydrothermal sedimentation.
The U/Th ratio can also be a proxy of redox conditions in paleo-seawater 34 .The U/Th ratio further confirms the impacts of hydrothermal activities, especially for phosphorite from the Gezhongwu Formation and dolomite from the Dengying Formation.Usually, normal sedimentary rock has U/Th < 1, while hydrothermal sedimentary rock has U/Th > 1 35 .The U/Th ratio (as shown in Table 4) in this study increases in the following sequence: carbonaceous shale from the Niutitang Formation (average value of 0.99, with values varying from 0.26 to 2.02) → phosphorite from the Doushantuo Formation (average value of 2.53, with values varying from 2.38 to 2.77) → dolomite from the Dengying Formation (12.26, with values varying from 5.33 to 16.45) → phosphorite from the Gezhongwu Formation (average of 18.23, with values varying from 5.99 to 56.81), showing the increasing impact from hydrothermal activities.The relationship between U and Th (Fig. 10) shows that the phosphorite from the Gezhongwu Formation and dolomite from the Dengying Formation are located in an area indicating paleohydrothermal sedimentation.
Evidence from rare earth elements Generally, positive Eu anomalies are widely used to indicate the involvement of hydrothermal materials.In this study, δEu values vary from high to low in the sequence of from the Dengying Formation dolomite (1.20, 0.97-1.57)> Gezhongwu Formation phosphorite (1.19, 1.03-1.44)> Niutitang Formation carbonaceous shale (0.91, 0.81-0.99)> Doushantuo Formation phosphorite (0.88, 0.83-0.98).Particularly, with the increase in the REE content in the bulk sediment of phosphorite from the Gezhongwu Formation, the enrichment of REEs is found to be closely related to the negative Eu anomaly (Fig. 11).Moreover, phosphorite from the Doushantuo Figure 9.The ratios of Fe/Ti, (Fe + Mn)/Ti and Al/(Al + Fe + Mn) of samples from the Zhijin and Weng'an phosphorite deposits.
Formation of the Weng'an deposit exhibits obvious negative Eu anomalies, indicating that REE-rich phosphorite from the Gezhongwu Formation developed with the involvement of hydrothermal fluid.The impact of hydrothermal activity continued until the deposition of the Dengying Formation.
A positive Y anomaly is also a geochemical indicator of hydrothermal sedimentation 36 .The average values of δY in this study vary from high to low from the Gezhongwu Formation phosphorite (2.22) > Dengying Formation dolomite (1.83) > Doushantuo Formation phosphorite (1.18) > Niutitang Formation carbonaceous shale (1.08); this sequence reflects that the phosphorite from the Gezhongwu Formation was formed with strong hydrothermal participation.Its highest REE content is then deduced to be closely related to hydrothermal activity.In the La/ Yb REE plot, Fig. 12 shows that most phosphorite from the Gezhongwu Formation is located in the intersection  environment: > 10) 39 .All of the samples in this study have (MgO/Al 2 O 3 ) × 100 > 10.However, the study area may experience changes between shallow sea and bathyal marine environments (Fig. 14).Generally, the seawater depth in the coastal zone is closely related to tectonic activity and sea level changes caused by transgression or regression 40 .The marginal shallow sea has lower salinity for the acceptance of freshwater from continental rivers.Moreover, the salinity in the near sea surface increases with the evaporation of seawater; therefore, the salinity of seawater is an important agent indicating variation in the seawater depth 41 .In addition, a deeper seawater depth means weaker dynamic conditions, and vice versa.In this sense, paleosalinity can also be used to reveal the paleoceanographic hydrodynamic environment.
The paleosalinity index of CaO/(Fe + CaO) (CaO/(Fe + CaO) < 0.2) represents low salinity, the ratio of 0.2-0.5 represents medium salinity, and high salinity is determined when CaO/(Fe + CaO) > 0.5 42 .The usage of the above index indicates that paleosalinity varied in different geological periods.The older strata, including phosphorite from the Doushantuo Formation and dolomite from the Dengying Formation, show high paleosalinity, with CaO/(Fe + CaO) > 0.5.The phosphorite from the Gezhongwu Formation is characterized by fluctuations in CaO/ (Fe + CaO) varying from 0.16 to 0.90, with an average ratio of 0.70, which indicates that the paleosalinity in this area varies widely.The CaO/(Fe + CaO) ratio of carbonaceous material from the Niutitang Formation varies from 0.01 to 0.23, indicating low and medium paleosalinities.The ancient seawater salinity in the Guizhou phosphorite area was generally high from the late Sinian to the early Cambrian, while from the Gezhongwu Formation to the Niutitang Formation, the salinity experienced a rapid decline, which may reflect the rapid rise in sea level, the deepening of sea water and other changes in the marine environment at that time.This may indicate the occurrence of special paleoceanographic events, such as violent transgression and regression or upwelling.The black shale unit at the top of the Doushantuo Formation is traditionally interpreted as transgressive deposits formed during sea level rise 43,44 .
The paleoceanographical hydrodynamics can also be revealed with element ratios such as Zr/Rb 45 .The grain size fractionating characteristics may be archived in the Zr/Rb ratio.The grain size sorting process during transport dominates the Zr and Rb abundance because Zr is preferentially enriched in the coarser grain size fraction, while Rb tends to be enriched in the finer grain size fraction.Thus, the Zr/Rb ratio increases with the strength of seawater energy 46 .In this study, dolomite from the Dengying Formation has Zr/Rb values from 0.84 to 10.00, with an average of 3.27, indicating a weak hydrodynamic environment.Zr/Rb in the phosphorite from the Gezhongwu Formation fluctuates greatly, varying from 0.65 to 65.45, with an average of 20.19, indicating a turbulent marine environment.However, both the carbonaceous shale from the Niutitang Formation and phosphorite from the Doushantuo Formation have Zr/Rb ratios of 3.29, varying from 0.88 to 7.49, and 1.34, varying from 0.42 to 2.80, respectively.Zr/Rb ratios indicate that the study area shows generally weak marine hydrodynamics from the late Sinian to early Cambrian, except during the Gezhongwu Formation, which showed a turbulent hydrodynamic environment and then weakened gradually to the Niutitang Formation.The recovery of paleoceanographic hydrodynamics is consistent with the results of petrological analysis of phosphorite from the Gezhongwu Formation.The phosphorite from the Gezhongwu Formation of the Zhijin deposit generally contains parallel bedding and cross-bedding, and some phosphorite contains wavy bedding, indicating that the phosphorite of the Zhijin deposit is subject to stronger hydrodynamic action during sedimentation.At the same time, the particle size of the collophanite in the phosphorite from the Gezhongwu Formation of the Zhijin deposit is larger than that from the Doushantuo Formation of the Weng'an deposit, indicating that the former developed higher textural maturity, which also indicates that the Zhijin phosphorite is subject to stronger mechanical action by seawater under shallower seawater depths than the phosphorite of the Weng'an deposit.

Tectonic setting analysis
The Zhijin phosphorite deposit in Guizhou is located at the southwestern end of the Central Guizhou Uplift in the Yangtze Craton.Most of the faults are linear ladder-like normal faults with structural characteristics similar to those in the ore-rich areas in the margin of the craton.The paleogeographic environment was a shallow sea with weak oxidation 47 .The Weng'an phosphorite deposit in Guizhou is located in the southeastern margin of the Yangtze Craton.The phosphorite of the early Ediacaran Doushantuo Formation was formed in a deep-water terrigenous marine environment 48 .In addition, the contribution of terrigenous weathering materials to the REEs in the phosphorites increased gradually from the early Ediacaran Doushantuo Formation to the early Cambrian Meishucun Formation 48 .Based on the water depth, the Zhijin phosphorite was subjected to late diagenesis and hydrodynamic actions over a longer time period, which promoted the adsorption of REEs by apatite from ambient seawater or pore water 49 .In contrast, the REE contents of Weng'an phosphorite deposits are much lower.Since the Weng'an phosphorite is part of the early Ediacaran Doushantuo Formation, the seawater was relatively anoxic with less biological activity.It lacked a sufficient supply of REEs and an appropriate depositional environment for accumulation.By comparison, the weakly oxidizing environment in which the Zhijin phosphorite deposit formed might have been an important factor controlling REE enrichment.
Furthermore, the REE-rich Zhijin phosphorite deposit may have been related to volcanic hydrothermal activity.Early Cambrian deep-sea hydrothermal vents have been identified in Zunyi, Guizhou.The major and trace element, REE, and sulfur isotopic compositions all indicate that a hydrothermal eruption occurred in this area in the early Cambrian 50 .These hydrothermal materials migrated to the shallow-water areas in Zhijin through upwelling, which may have caused the REE enrichment of the Zhijin phosphorite.However, the Weng'an phosphorite was formed in the late Ediacaran, i.e., earlier than the hydrothermal vent activity.These may be the reasons for the low REE contents of the Weng'an phosphorites.Oxygen isotope geothermometry also indicates that the REE enrichment of the Zhijin phosphorite is related to higher-temperature hydrothermal fluids 48 .
In summary, the tectonic activity on the craton margin provides a suitable depositional environment and material supply for the formation of REE-rich phosphorite.During tectonic activity, the input of volcanic

Figure 7 .
Figure 7. PAAS-normalized rare earth element distribution patterns of samples from the Zhijin and Weng'an phosphorite deposits.(a) Phosphorites from the Doushantuo Formation.(b) Dolomite from the Dengying Formation.(c) Phosphorites from the Gezhongwu Formation.d Carbonaceous shale from the Niutitang Formation.

Figure 8 .
Figure 8. ε Nd (t) − ( 87 Sr/ 86 Sr) i diagram of samples from the Zhijin and Weng'an phosphorite deposits.The upper crust data are from Jahn et al.28 .The ancient lithospheric mantle data are from Ma et al.29 .

Figure 10 .
Figure 10.Co-Ni-Zn and U-Th diagram of samples from the Zhijin and Weng'an phosphorite deposits.

Characteristics of major elements The
K 2 O + Na 2 O) contents are very low.TFe 2 O 3 values vary in the range of 0.12-1.34%,withan average value of 0.59%.In addition, the LOI values in samples from the dolomites of the Dengying Formation are very high, with an average value of 32.33% and values varying from 24.32 to 45.41.The phosphorites from the Gezhongwu Formation contain very high SiO 2 , with an average value of 44.97% and values varying in a wide range from 24.04 to very 65.90%.These phosphorites contain high TFe 2 O 3 , with an average value of 6.14%.The CaO content is intermediate in comparison with rocks from other formations, with an average value of 20.68%.However, MgO (average value of 1.50%), Al 2 O 3 , TiO 2 and (K 2 O + Na 2 O) are low.A relatively high LOI is observed in phosphorites from the Gezhongwu Formation, with an average value of 13.59% and values ranging from 10.52 to 18.61%.

Table 2 .
The original data of major elements (wt%), trace elements (ppm) by ICP-MS and rare earth elements (ppm) by ICP-AES in this study.LOI loss on ignition.TFe 2 O 3 = FeO/0.8998.

ppm) 147 Sm/ 144 Nd 143 Nd/ 144
The carbonaceous shale from the Niutitang Formation of the Zhijin deposit is characterized by high contents of SiO 2 and Al 2 O 3 , with average values of 63.16% and 17.42%, respectively, with values varying from 60.84.25 to 65.47% and 15.31 to 19.20%, respectively.TFe 2 O 3 is higher, with an average value of 4.86% and values ranging from 2.72 to 6.40%.Obvious scattered spots of pyrite are observed on the rock surface.However, there are low contents of CaO and MgO in these shales, with average values of 0.58% and 2.16%, respectively; however, the (K 2 O + Na 2 O) content is high, with an average value of 4.42%.

Table 4 .
The statistics of the major elements, trace elements and rare earth elements compositions in this study.